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Abstract 

What is the connection between how the dark matter was produced in the early 
universe and how we can detect it today? Where does the WIMP miracle come 
from, and is it really a "WIMP" miracle? What brackets the mass range for thermal 
relics? Where does {av) come from, and what does it mean? What is the differ- 
ence between chemical and kinetic decoupling? Why do some people think that 
dark matter cannot be lighter than 40 GeV? Why is bh such a popular annihilation 
final state? Why is antimatter a good way to look for dark matter? Why should 
the cosmic-ray positron fraction decline with energy, and why does it not? How 
does one calculate the flux of neutrinos from dark matter annihilation in a celes- 
tial body, and when is that flux independent of the dark matter pair-annihilation 
rate? How does dark matter produce photons? — Read these lecture notes, do the 
suggested 10 exercises, and you will find answers to all of these questions (and 
to many more on what You Always Wanted to Know About Dark Matter But Were 
Afraid to Ask). 



* Lecture Notes for TASI 2012: Theoretical Advanced Study Institute in Elementary Particle Physics 
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Introduction 



These are lecture notes for a series of four lectures on astrophysical probes of dark 
matter I delivered in June 2012 at the Theoretical Advanced Study Institute in Ele- 
mentary Particle Physics (TASI) Summer School at the University of Colorado, Boul- 
der. The aim of this set of lectures is neither to present in detail particle dark matter 
models nor to focus on the technical aspects of experiments or the (alas, necessary!) 
understanding of astrophysical backgrounds; rather, I try to convey and to work out 
order-of-magnitude estimates that can be applied to a variety of particle dark matter 
models and physical situations. Practice with these estimates might be perhaps appre- 
ciated by the theory-minded scholar who is not necessarily keen on the fine print of 
experimental setups or on the complications of astrophysical processes. 

The first lecture is devoted to an introduction to the cosmology of particle dark 
matter, specifically in connection with the indirect detection of dark matter: what are 
the relevant/expected energy scales, particle products and rates? The second lecture 
discusses the thermal (or non-thermal) processes associated with the production of 
dark matter in the early universe, a few exceptions to the "standard lore", and the 
process of kinetic decoupling. The third lecture introduces indirect detection of dark 
matter, and details on charged cosmic rays produced by dark matter in the Galaxy. 
Finally, the fourth lecture discusses gamma rays and high-energy neutrinos from dark 
matter. Ten simple exercises are scattered through the lectures. The interested Reader 
might find them useful practice to master the material being discussed. 

The reference list is vastly incomplete, and it primarily egotistically contains self- 
citations: I consider it part of a good scientific education to go the extra mile and 
find relevant papers on given topics of interest. This will therefore be Exercise #11! 
Enjoy! 
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Lecture 1: Particle Dark Matter: zeroth-order Lessons 
from Cosmology 

The fundamental (elementary) particle nature of dark matter can be probed with 
the detection of photons, neutrinos or charged cosmic rays that are produced by, or 
affected by, dark matter as an elementary particle. The key processes are: 

(a) the pair annihilation of dark matter particles (which we shall generically indicate 
here with the symbol x), producing Standard Model (SM) particles in the final 
state: x + X^ SM; 

(b) the decay of dark matter particles into SM particles: x SM; 

(c) the elastic scattering of dark matter particles off of SM particles: x + SM 
X + SM. 

Other processes might exist, but are less common in the literature and in model build- 
ing, and we will not entertain them here. It is important to note that none of the 
processes (a-c) listed above is bound to necessarily occur in any particle dark matter 
model: for example (b) does not occur if the dark matter is absolutely stable, and 
(a) and (c) can be highly suppressed, or even not occur at all, if the coupling of the 
dark matter sector to the SM is suppressed, or if the dark matter sector is somehow 
"secluded", or if the dark matter is not its own antiparticle, and there is no anti-dark 
matter around. 

There exist, however, reasons to be optimistic with respect to the prospect of de- 
tecting non-gravitational signatures from dark matter: firstly, some of the best mo- 
tivated (from a theoretical standpoint) extensions to the SM encompassing a dark 
matter candidate x predict a coupling of x to SM particles that would entail processes 
(a)-(c) or (b) at some level; secondly, there exist "phenomenological" reasons, chiefly 
the so-called WIMP miracle (to be reviewed in what follows), that imply the occur- 
rence of some of the processes above for models where the observed abundance of 
dark matter is connected (via thermodynamics and cosmology) to its particle nature. 

There are three key ingredients to understand indirect dark matter detection at a 
qualitative level, and to be able to make quantitative predictions: 

1. Production rates of the relevant SM particles ("messengers"); this is related to 
the pair-annihilation or decay rate of the dark matter particle; 

2. Energy scale of the SM messengers: this is set by the mass of the dark matter 
particle (or by its momentum, for processes of t5^e (c) above); 

3. Annihilation products: this largely model-dependent ingredient details on which 
SM particles are produced by the dark matter particle. 
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The rate re±,p,^,iy,... for a given SM messenger (i.e. the flux of such particle species 
per unit time from a unit volume V containing dark matter particles) is generically 
the product of three factors: (1) the number of dark matter particle pairs (or of dark 
matter particles, for decaying %) in the volume V times (2) the pair annihilation (re- 
spectively the decay) rate, times (3) the flux of SM particles per annihilation (decay) 
event. In formulae: 

TSM, ann = ^^V^ >< (^^) >< (^SM, ann) , 

TsM, dec = (y"~^^) ^ (-7^) ^ 

Interestingly for the present discussion, many key quantities (pdm, "^x' • - ^ 
potentially connected to how the dark matter was produced in the early universe. The 
dark matter production mechanism in the very early universe is therefore a great start- 
ing point both for model building and for eye-balling the relevant indirect detection 
techniques and for setting constraints. 

The one quantity from cosmology which is important to have in mind is the average 
dark matter density in the universe, 

pDM ^ 0.23 ■ Pcrit = ^DM ° 



8nG 



N 

It is useful in many social (as well as anti-social) situations to have on the tip of 
your tongue the value of this latter quantity both in "astronomical"^ and in "particle 
physics" units: 

Mpc cm'^ 

From the "astronomical" units, we learn for example that clusters of galaxies, the 
largest bound dark matter structures in the universe, have typical over-densities^ of 
10^, since they approximately host hundreds to thousands of galaxies, whose mass is in 
the ~ 10^^ Mq range; from the "particle physics" units, we learn that in our particular 
location in the Milky Way, where pdm ~ 0.3 GeV/cm^, the over-density is a factor of a 
few larger than in a typical cluster 

A successful framework for the origin of species in the early universe is the paradigm 
of thermal decoupling (see e.g. [1] and [2]). This framework encompasses for exam- 
ple the successful predictions of recombination and of Big Bang nucleosynthesis, and 
it describes in detail the process of cosmological neutrino decoupling. In short, ther- 
mal decoupling consists of the relevant particle interaction rate F, initially, at high 
temperatures, much larger than the Hubble expansion rate H, falling to a "freeze-out" 



^Particle physicist: always a good idea to talk to astronomers; for example, I found my wife that 
way! 

^An "over-density" is a region whose average density is larger, by a certain (over-density) factor, 
than the overall average density. 
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point where F ~ iJ; after this point in time/temperature, the particle species simply 
"redshifts" its momentum and number density. 

In natural units, we can think oiT = n ■ a, with n a particle number density and a 
an interaction cross section. As statistical mechanics kindly teaches us, the equilibrium 
number density of a particle of mass m in a thermal bath of temperature T has two 
asymptotic regimes: 

rirei ~ for m<r, 

^non-rei ~ (mTf^'^expf-^) for m > T. 



The right-hand side of F ~ if, i.e. H{T), comes from general relativity, and specifically 
from Friedmann's equation: 

2 SttGn 
H = —p. 

In the radiation dominated epoch (i.e. T > 1 eV), 

with g the number of relativistic degrees of freedom ig = 2 for photons). To a decent 

degree of approximation, and recalling that Mp = I/^/S-kGn, in order to eyeball when 
thermal decoupling occurs you can take H ~ T"^ /Mp. 

Let's put all of this in practice, and estimate the temperature of neutrino freeze- 
out. We estimate the relevant scattering cross section in the Fermi four-fermion 
contact interaction approximation, and we take E ~ T^, so that a ~ C^T^ (where 
Gp ^ 10~^ GeV~^ is Fermi's constant). We thus have 

n-a^H ^ T^GlT^ = T^/Mp. 

We thus have 

n = (GlMp)-^/^ ~ (10-^° X 10^^)-^/^ GeV ~ 1 MeV. 

Among various things to be happy about, we cheerfully verify that the T,y we found 
indeed satisfies T„ » m^, which we have implicitly assumed for the form of n{T), and 
we learn that neutrinos are hot relics, not because they are particularly attractive, but 
because they freeze out while they are relativistic. 

Now let's calculate the relic density for a dark matter particle x- L^t me introduce 
the notation Y = n/s where n is a number density and s is the entropy density. In an 
iso-entropic universe, s • = constant, where a is the universe's scale factor. Y ~ na^ 
is thus a "comoving" number density. If no entropy is produced, ytoday = freeze-out- In 
the case of hot relics, like SM neutrinos, 

-'freeze— out 
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and 

^today "^today ^ freeze— out- 

For example, for the SM neutrinos, the fraction of the universe's critical density times 
h'^ (where h is today's Hubble constant in units of 100 km/s/Mpc - in practice h'^ ~ 
0.5) in a SM neutrino species is 



Pcrit 91.5 eV 

While the normalization depends on the relevant cross section, it is a general fact that 
a hot relic's thermal relic abundance scales linearly with the relic's mass. For a weakly 
interacting dark matter particle, requiring that the thermal dark matter density be less 
or equal than the observed matter density leads to the so-called Coswik-McClelland 
limit [3] on the mass of a hot dark matter relic. 



Exercise #1: Calculate Tf.o. for the pp annihilation reaction (you can 
use a ~ m;^^) and estimate the relic proton/antiproton density; is 
this a hot relic problem? Compare what you find with the observed 
"baryon asymmetry". 



A cold relic is one for which the freeze-out temperature is much lower than the mass 
of the particle, which thus decouples in the non-relativistic regime. An illustrative 
example of a cold relic is a "heavy" neutrino, with a mass tun 1 MeV. In this 
case, the appropriate asymptotic form for the equilibrium number density is the non- 
relativistic limit 

Again, the condition na H yields 

Mp ■ a 

Let us call m^/T = x; when dealing with cold relics, we are thus working in the x 3> 1 
regime. We can re-cast the condition n ■ cr ~ i7 as 



2,3/2 x^-Mp-a' 



We thus need to solve 



X ■ e'"" = 10"^^ (2) 

m^-Mp-a 102 ■ 1018 ■ 10-6 ' ^ ^ 

where I've substituted for the nominal values of an "electro-weak interacting" relic, 
with a ~ Gprn"^ and ~ 10^ GeV. Numerically, for the range 10"i°...10-2o for the 
right-hand side of Eq. (2), the resulting a^f.o. — 20. ..50. Now, 

_ • %(T = To) _ np 

He Pc -1-0 
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with To = 2.75 K ~ 10 eV. Since for an iso-entropic universe aT ~const, 

To' Tt. 

we have 

where I used Eq. (1) in the last step. The equation above can be then cast, plugging 
in the numbers for the various constants, as 



xi,o. flQ-^ GeV"^ 



0.2/ 20 V a 



(3) 



a relation that many refer to as "miraculous". Often, Eq. (3) is quoted with the 
thermally-averaged product of the cross section times velocity {av) (we will under- 
stand why this is, and what a thermal average is, in the next lecture), instead of the 
simple cross section a. Since v ~ c/3 for x ~ 20, one has 



cm Oft cm^ 



(av 



) ~ 10-^ (3 X 10-28 GeV^ cm^) 10^° — = 3 x 10-^^ 



Exercise #2: Convince yourself that v ~ c/3 for x ~ 20. 



The {av) ~ 3 x 10 cm^/s is a "magic" number definitely worth keeping in mind! 

Is the "magic number" we just found unique and peculiar to the electroweak scale? 
Not at all! Let us remind ourselves which ingredients we used to get the "right" relic 
density: 

(i) the condition for having a cold relic, ■ a ■ Mp 3> 1; 

(ii) a cross section a ~ 10"*^ GeV^^. 
Now, suppose that the cross section be 

a 

mi 

A. 

with g some coupling. Condition (ii) essentially enforces that 

2 

g ~ 



10 TeV 

independently of which scale is at. Now go back to condition (i) : this reads 

1 ml 10^ GeV^ 



Mp > 



m^ ■ a ■ {g"^)^ m 



X 
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if condition (ii) holds. Therefore, condition (i) impUes, with condition (ii), that 
rriy^ ^0.1 eV. Therefore, thermal freeze-out giving the "right" relic abundance is not 
peculiar to the electroweak scale, as reiterated recently in the literature (see e.g. the 
"WIMPless" miracle of Ref. [4]). However, since 

^Ew - GlTl, ~ ^^(^) ~ 10"' GeV-^ 

the electroweak scale is quite a "natural" place (whatever natural means) for the mir- 
acle to occur! 

Is there any upper limit to the particle dark matter mass in the cold thermal relic 
scheme? Indeed there is! The coupling constant g cannot be arbitrarily large (a 
condition that can also be rephrased in terms of a unitarity limit in the partial wave 
expansion [5] ; note that caveats to the unitarity argument do exist, and this limit can 
be evaded! (I suggest you read Ref. [5] and think about how to do that)). Roughly, 



ml 



which implies 

— > 10"'' ^ 

0.2 ~ 

Therefore, demanding < 0.2 implies 



0.2 ~ 47r 



m 



X 



120 TeV 



2 



<1, 



or < 120 TeV. 



Is there, similarly, a lower limit in the cold thermal relic scheme? We commented 
above on the general limit, for arbitrarily low cross sections, ^0.1 eV. But sup- 
pose now we have in mind a particle that interacts via electroweak interactions, for 
example, again, a massive neutrino with a ~ Gj^ m^. In this case 



X 

2 



2 10-*^ 1 /10GeV\- 
2^0 1 ■ ~ 1 



This implies that m^ > 10 GeV for weakly interacting massive particles (WIMPs) - a 
limit known in the literature as the Lee-Weinberg limit [6] . 

Fig. 1, from Ref. [7], illustrates the thermal relic density of a weakly interacting 
massive particle as a function of the particle's mass. The cross section is assumed to 
be of the form 



2 



cr ~ (4) 

[s — m^/j + m^, 

with s the total center of mass energy squared. The mass of the mediator Z' is taken 
to be 10 GeV, 91.2 GeV (the Z mass) and 1 TeV. The asymptotic hot and cold relic 
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log (myCeV) 



Figure 1 : The thermal reUc density of a reUc that pair-annihilates with the cross section 
of Eq. (4), for three values of m^/. From Ref. [7]. 

behaviors are clearly visible and match the predictions we made above: a hot relic 
density scales linearly with mass, a cold relic with ^ mz' has i7 ~ 1/cr ~ m^, and 
a cold relic in the regime where ^mz' has Vt ~ m^i/rn^^. 

In theories with a more complicated dark matter pair-annihilation cross section 
than what appears in Eq. (4), for example for the lightest neutralino in the minimal 
supersymmetric extension of the Standard Model (MSSM), various effects blur the 
simple connection between mass and cross section/relic density, resulting in a much 
larger spread of results. The general feature of a lower limit of about 1-10 GeV is, 
however, rather resilient, as shown in Fig. 2, from Ref. [8], where I scanned generously 
over the relevant MSSM parameter space. The black dots indicate points with tan /3 
(the ratio of the vacuum expectation values of the two Higgses in the MSSM) fixed to 
50, while the red points have tan /3 = 5. I invite the SUSY-afecionado to read my paper 
for details of the scan. The x-axis indicates the lightest neutralino mass, while the y- 
axis is the thermal relic density calculated in a standard cosmology without entrooy 
injection or a modified Hubble expansion rate. I stretched the relevant parameters as 
hard as I could to obtain the "optimistic limit" green line. That line shows that for 
Vlh"^ ~ 0.1 the lightest neutralinos I find are in the few GeV range. 

Let us now ask the general question: how light can a WIMP be? So far, we assumed 
an iso-entropic universe. Suppose at some point after a WIMP has frozen out and is 
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Figure 2: The thermal reUc density of the Hghtest neutraUno in the MSSM, as a func- 
tion of the neutraUno mass. From Ref. [8]. 



thus decoupled from the universe's thermal bath, the entropy density changes from 
s — )■ 7 ■ s, 7 > 1 from e.g. decaying relics (such as relic gravitinos, moduli, . . . ) or 
from a first order phase transition (see e.g. Ref. [9]). Then 

^today ^^^^ and n^^^. 

7 7 

For a sufficiently large 7, the relic abundance of almost any over-abundant relic WIMP 
can be "diluted" enough to match the observed dark matter density. For example, in 
supersymmetry the lightest neutraUno can be almost arbitrarily light as long as it is 
bino-like (to prevent an excessively light associated chargino) and if sfermions are 
sufficiently heavy to suppress energy loss mechanisms in stars (for example e+e 
Xx)- The additional requirement is of course that the entropy injection happen at a 
temperature smaller than the neutraUno freeze-out, which sets a (weak) constraint 
on the neutraUno mass [8] : in order to maintain the successful predictions of light 
elemental abundances, entropy injection cannot happen too close to the era of Big 
Bang Nucleosynthesis. 
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Lecture 2: WIMP Relic Density, a Closer Look 

The dark matter literature is flooded with the symbol {<jv), short for the zero-temperature 
thermally-averaged pair-annihilation cross section times velocity - but do we really un- 
derstand what this symbol indicates, and where it comes from? How is the "thermal 
average" defined? What is v7 is it a relative velocity? but a relative velocity is not 
Lorentz-invariant! etc... 

The starting point for a closer look at the WIMP relic density is the Boltzmann 
equation, that can be symbolically cast as [1]: 

L[f] = C[f], (5) 

where / = f{p,x,t) is the phase space density, L is the Liouville operator describing 
the change in time of the phase space density, and C is the collision operator describing 
the number of particles per phase-space volume lost or gained per unit time. For those 
(like me) who need a refresher on the Liouville operator, its non-relativistic form reads 

while in its covariant form it reads 

(9 (9 

In a homogeneous and isotropic cosmology (also known in the trade's slang as a 
Friedman- Robertson-Walker universe) , 

/(|pl,t) or, equivalently : f{E,t). 

Also, L simplifies to 

. df (2 2 9f 

^f^] = ^9^-a"^^ dE- 
We are interested in particle number densities, defined by 



spin 



We will thus take Eq. (5) and consider (calling g the number of spin degrees of free- 
dom) 

where we introduced H = a/a and integrated by parts using 

1 d / o X dn 
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Cleaning up the right-hand side of the Boltzmann equation is a bit messier and I 
recommend the classic paper by Gondolo and Gelmini, Ref. [10]. For definiteness, 
let us consider a process of the type 1 + 2—^3 + 4 where we are interested in the 
number density of species 1, and where we assume that species 3 and 4 are in thermal 
equilibrium. The right-hand side of Eq. (5) can then be cast as: 

= -{(r- VM0i) {riiUi - n^X*^) , 

where ni,2 are the number densities, while nl\ indicate the equilibrium number den- 
sities, and where 

f 

indicates the invariant, unpolarized total cross section for processes 1 + 2 — >■ any final 
state / in thermal equilibrium, and where, finally, the "M0ller velocity-^" is defined in 
the following covariant form: 



El E2 

A couple of comments: 

1. Notice that vm0\ ^1 "^2 is a Lorentz invariant quantity; 

2. Notice that in the rest frame of 1 (or 2; what we can think of as the "lab frame"), 
VM01 Vrei = \vi - V2\, where e.g. vi = pi/Ei etc. 

Last ingredient: the thermal average: this is defined by the expression: 



Exercise #3: Evaluate the denominator of Eq. (6) for nii = 1712. 



The diligent Reader who carried out the exercise proposed above found that the de- 
nominator of Eq. (6), for rrii = 1112 = m (the pair-annihilation case relevant for us) 
reads 

I e-^^l^e-^-l^ d^vi dh^ = {47rm'TK2 (^))' , 

where K2 is the modified Bessel function of the second order The numerator reads, 
instead, 

J <y ■ ^M.ie-^^/^e-^^/^ d'pi d'p2 = j^^ <y{s - Am')^sKi d., (7) 



^Usually the M0ller velocity is defined as WM01 — {{vi — ^2)^ — {vi x V2Y 
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a "convolution" of the cross section evaluated at s — 4m^, where s is the center-of- 
mass total energy, and a temperature-dependent thermal kernel. This form is key to 
understand important caveats to, e.g. the magic relation that implies that {av) = 
3 X 10"^^ cm^/s gives the correct thermal relic density. From now on, I will suppress 
the subscript M0I and intend always that v — )■ vm0\- 

Caveats to the Standard Story 

A classic paper that discusses departures from the vanilla calculation of the WIMP 
relic density described in the previous lecture is Ref. [11] by Griest and Seckel, "Three 
exceptions in the calculation of relic abundances". The three memorable exceptions 
are: 

1. Resonances; 

2. Thresholds; 

3. Co-annihilations. 

Resonant annihilation through a particle with the right quantum numbers and a mass 
niA — 2m^, found for example in the so-called "funnel" region of the (soon to be gone, 
thanks LHC!) minimal supergravity/constrained MSSM model, can be relevant either 
if > or if < In the first case, the cross section peaks at s = 

and is thus most relevant at temperatures Trcs — rn\/ (Qm^). 

Exercise #4: Show that (s) ~ 4m^ + 6m^T. 



If Tres — Tf.o. the resonance is extremely important at freeze-out, and hence for the 
thermal relic density; the pair annihilation cross section today, the one that we care 
about for indirect dark matter detection rates, will however be potentially much lower! 
In this case, the freeze-out cross section might be close to {av) = 3 x 10"^^ cm^/s, but 
(potentially) the T = cross section {av)o ^ 3 x 10"^^ cm^/s. 

If < ■mA/2, in the integral of Eq. (7) the cross section is always maximal for 
T = 0, the resonance can be (or not) subdominant at freeze-out, but we are in the 
(lucky if one wants a signal, unlucky if one wants to hide it!) circumstance where the 
T = cross section(crf )o 3> 3 x 10~^^ cm^/s. It must be noted that this discussion 
has model-dependent caveats: for example, in supersymmetry the lightest neutralinos 
are Majorana particles, and in a purely s-wave annihilation (at T = 0) a pair of 
neutralinos are in a CP-odd state. Therefore, for example, the pair annihilation via a 
CP-even particle (such as a neutral CP-even Higgs) cannot contribute to the T = 
pair annihilation cross section! 

Thresholds affect the relation between the freeze-out and T = pair-annihilation 
cross sections in an obvious way: the cross section in Eq. (7) suddenly increases as, 
e.g. s > Amf, where mt is the particle in which pair our annihilating particle can go. 
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XX ~^ tt (think e.g. of < rrit where t is the Standard Model top quark). Thresholds 
therefore always imply (erf )f.o. > {pv)^. 

Co-annihilation occur for particles whose freeze-out process is tangled with that 
of other particle species with a close enough mass so that the two freeze-out episodes 
are inter-connected. A necessary condition is that this second co-annihilating species 
2 have a mass such that at freeze-out the Boltzmann suppression of its equilibrium 
number density is not dramatic. In formulae, we want m2 — mi < Tf o.. In this case 
the relevant cross section is an "effective" cross section that includes the appropriately 
Boltzmann-weighed contribution from (A^) co-annihilating particles, i.e., with obvious 
notation (if this is not obvious see the extensive and clear discussion of Ref. [12]) 



In the equation above, Amj indicates the difference in mass between particle i and the 
lightest particle (to which i, eventually, decays) . Note that the denominator counts the 
effect of the additional degrees of freedom, suitably weighed. 

Coannihilation comes in two varieties, that I like to call "parasitic" and "symbi- 
otic". If the additional degrees of freedom annihilate "less efficiently" than the parti- 
cle whose number density we are interested in, then the coannihilating particles will 
have a parasitic effect, and produce a smaller effective pair-annihilation cross sec- 
tion. This is the typical case with, for example, the lightest Kaluza-Klein excitation of 
universal extra dimensions (UED; for a review see my own review! Ref. [13]). UED 
has a very compressed spectrum of particles (Kaluza-Klein modes of Standard Model 
particles, whose mass differs from the compactification scale by loop corrections or 
corrections of the order of the Standard Model particle mass) above the stable dark 
matter particle candidate. This large collection of particles brings many additional 
effective degrees of freedom which outweigh the corresponding additional contribu- 
tion to the pair-annihilation cross section, rendering UED a prototypical example of 
parasitic coannihilation. 

An example of "symbiotic" co-annihilation is provided by a nimble particle that an- 
nihilates efficiently and that doesn't carry a large number of degrees of freedom. An 
example is co-annihilation of the lightest neutralino of the MSSM with the scalar part- 
ner of the tau lepton, or "stau" (but other equally good examples are co-annihilation 
with charginos or stops or other sfermions) . Unfortunately, collider searches indicate 
that the days of the "stau coannihilation region" may be counted [14]. 

So far we have dealt with exceptions to the left-hand side of the Boltzmann equa- 
tion (or of its short-hand version T = n ■ a = H). What happens if we fiddle around 
with the right-hand side instead, i.e. with H, the expansion history of the universe? To 
avoid self-promoting my own papers again, I will point you to the following example: 
a cosmology with a "quintessence" field that provides a dynamical dark energy term, 
whose impact on the relic density of WIMPs was first studied by Salati in Ref. [15]"*. 

'^See also Ref. [16]. Couldn't resist. 
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Let (j) be the quintessence field, a spatially homogeneous real, scalar field. The field 
energy density and pressure are 

An example of a suitable potential that exhibits the desired "tracking" behavior (for 
appropriate initial conditions), i.e. whose energy density tracks dynamically the dom- 
inant energy density component, is 



The field's equation of state w = P4,/ moves from w = — 1 in the "kination" phase, 
where the kinetic energy term dominates, to w = +1 in the "cosmological constant" 
phase, where V dominates. Tracking helps explain the coincidence problem = 
Vt^ ~ VIm, although fine-tuning is not eliminated (as it creeps back in via the field's 
initial conditions) . 

Noticing that ~ a~-^^^^'"\ in the kination phase p ~ and therefore the uni- 
verse is kination-dominated as sufficiently early times, with 

rp2 rj-i 

H , 

Mp Tkre 

where Tkre stands for the temperature of kination-radiation equality. To be relevant 
for the relic density of a particle species decoupling at T = Tf „ , kination must dom- 
inate before and at freeze-out, hence Tkre > Tf o.. However, to avoid disrupting Big 
Bang Nucleosynthesis, we must also require that Tkre < Trbn ~ 1 MeV. 

In a kination-dominated universe, freeze-out works, schematically, exactly as we 
described in the previous section, and 



T^3 



^^i^o. 

- Mp.p/'-\Tl^ 

but now the freeze-out condition reads 

T^ T 
ni.o.{(J v) ~ — . 

Mp 1 KRE 

We therefore have that 

rii.o. 1 Tf.o. 



Tl^ Mp{av)TKRE' 
To first order, the enhancement factor of the thermal relic density in the presence of 
quintessence above the standard thermal relic density is thus 

Q^tandard Tkre ~ 20 TbbN 100 GeV 
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With more accurate calculations the enhancement factor is found to be potentially as 
large as 10^ [15, 16]. 

The dark matter production mechanism can naturally be non-thermal, or it can 
arise from an asymmetry (as is probably the case for the origin of baryonic matter) . 
Suppose a particle species ip, with > is produced in the early universe with an 
abundance Q,^, and that ip decays to x, which is the stable dark matter particle, at a 
temperature when x is out of equilibrium. The relic density that x will then inherit (up 
to contribution from the decay of other particle species and from thermal production 
etc.) is simply 

where the ~ sign indicates that additional effects (such as some entropy production 
in the decay process) can enter. 

As you calculated in Ex. #1, the thermal abundance of protons and antiprotons 
is almost ten orders of magnitude smaller than the observed abundance. If an asym- 
metry is present, then the observed proton density can be inherited entirely from the 
asymmetry itself. The same could hold for the dark matter sector Many have envi- 
sioned the possibility of explaining the coincidence ^bm — (where DM is dark 
matter and B is baryonic matter) by postulating that perhaps n^u ~ ub and that 
~ 5 GeV, with some mechanism that couples an asymmetry in one sector to the 
other sector, or many variants on this theme. 

It is interesting to ask the question of whether any indirect dark matter detection 
signal could exist if dark matter originates from an asymmetry. In principle, if no 
anti-dark-matter is generated (and if primordial pair-annihilation got rid of all of it), 
then it's very hard to get any indirect signals; however, if dark matter and anti-dark 
matter oscillate into each other (for example because of a AA^^ = 2 operator, such as 
a mass term, say mMXX ) then oscillations can populate the anti-dark matter content 
of the universe and residual annihilation occur [17]. Unless a symmetry in the theory 
explicitly prohibits such operators, oscillations generically occur, since after all we only 
need 



Of course many other interesting and witty dark matter production mechanisms 
exist: I recommend you explore at least the classic ones that pertain to sterile neutri- 
nos and to axions (see e.g. Ref. [1]). 

Kinetic Decoupling 

The dark matter thermal history in the very early universe is important not only for 
the calculation of the particle's relic density, but potentially also for the formation of 
matter structure in the universe, especially for (cold) WIMPs. In the early universe, 
elastic scattering processes such as xf ^ xi\ where / is a Standard Model fermion. 
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keep the dark matter particle in kinetic equilibrium even after chemical decoupling 
(i.e. when T-^^^x ^ H). The reason is that the target densities for the processes that 
keep the dark matter in chemical versus kinetic equilibrium are vastly different after 
chemical decoupling: 

XX^ ff r = rinon-rcl " 0" 

Xf ^ Xf r = r2,ei ■ o- 

with nnon-rei expoueutially suppressed! Let us now estimate the kinetic decoupling 
temperature of a WIMP. We shall assume a default electroweak cross section 

We have to account for the fact that the WIMP is non-relativistic after chemical decou- 
pling and that momentum transfer between the WIMP and the thermal bath becomes 
"inefficient", in a sense to be made quantitative with a simple estimate. The typi- 
cal momentum transfer per collision is 5p ~ T, while the WIMP momentum in the 
non-relativistic regime satisfies the relation 

^ T 

2m 

so that p ~ ^myT. Momentum transfer is a stochastic process, so it takes 



p J f^x^ ^x 

collisions to establish kinetic equilibrium. 

To calculate kinetic decoupling, we thus want to compare 

f rr.3 r^2rr2 T 



We thus find 



p J M 



I \^^^ / 777,^ 

k.d. 



Tvri ~ ( , ~ 30 MeV ' ^ 



Mp-GlJ VlOOGeV. 

What does this imply for structure formation? Roughly, the cutoff scale of the matter 
power spectrum will correspond to the size of the horizon at kinetic decoupling, so 



47r / 1 \^ _ , ,/10 MeV^ ^ 



Mcutoff ~ ^ 7777^ PDM(Tkd) ~ 30 M, 



More precisely, the cutoff scale is set by the largest of the free-streaming versus acous- 
tic damping scale, see e.g. [18] for a nice review. For typical WIMPs, we thus find 
that these "protohalos" (which correspond to the first structures that gravitationally 
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Figure 3: The correlation between the spin-dependent dark matter-proton cross sec- 
tion and the kinetic decouphng temperature (upper panels) and the small scale cut- 
off mass (lower panels) in the MSSM (left panels) and in UED (right panel). From 
Ref. [20]. 



collapse in the early universe) have a mass comparable to the Earth mass (i.e. roughly 
10~^ Mq). However, in specific theories the range of variations can be very significant 
[19]. This fact has potentially important consequences for indirect detection, as it 
feeds in the problem of calculating the boost factor (i.e. the additional contribution to 
the count of pairs in a halo on top of the smooth halo component) from substructure, 
and for the dark matter "small-scale" problem. 

Is there anyway to probe the size of the dark matter small-scale cutoff? In Ref. [20] 
we pointed out that if / = g (i.e. the particles the dark matter scatters off of are 
quarks), the processes relevant for kinetic decoupling are exactly the same as those 
participating in dark matter direct detection. Therefore, in principle, one could corre- 
late Mcutoff with cTdircct dct- There is, however, an important caveat: usually, T^d < Tqcd, 
the latter symbol indicating the temperature corresponding to the QCD confinement 
phase transition. After confinement, the number density of hadrons is negligible com- 
pared to light leptons, and the latter dominate and control kinetic decoupling. If, 
however, one has a dark matter theory where "quark-lepton universality" holds, then 
a correlation is expected, and indeed is found - both for UED and for supersymmetric 
dark matter [20] as illustrated in fig. 3. 
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Lecture 3: Indirect Dark Matter Detection 



What have we learned thus far about annihilation processes that we could use to 
detect non-gravitational signals from particle dark matter? 

- {av) ~ 3 X 10"^® cm^/s is a good "magic" benchmark; 

- the magic number above is not WIMP-specific, and is independent of mass, to 
first order; 

- there exist numerous caveats, both from the particle physics side (coannihilation, 
resonances, thresholds,...) and from the cosmology side (quintessence, non- 
thermal production, asymmetry...) to that magic number 

What about dark matter decay? If dark matter is unstable with a lifetime well in 
excess of the age of the universe, the decay products would also be a great way to 
detect non-gravitational signals from dark matter! From a theoretical standpoint, a 
GUT-scale or Planck-scale dark matter number-violating operators should be generic. 
For example for a dimension-5 operator, 



the resulting lifetime 



would potentially result in an impact on BBN, but it would imply an excessively short- 
lived dark matter candidate^. 

For a dimension-6 operator things look more interesting. 



with a lifetime 



which turns out to be a very interesting lifetime range to explain the Pamela positron 
excess and for searches for dark matter with gamma rays, as we shall see later on. 

For definiteness, let me however concentrate on dark matter annihilation. Let us 
now try to corner the key ingredients to make predictions for indirect searches for 
dark matter First and foremost: what do we know about the dark matter particle 
mass, which sets the energy scale for the particles produced in an annihilation event? 

^As a supplementary exercise, make sure you refresh your memory about "natural units" and how 
to convert energies into times, and convince yourself of the overall factor of 1 s in Eq. (10) 




(10) 
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We saw earlier that for WIMPs a useful lower limit is provided by the Lee-Weinberg 
bound at about 10 GeV, while unitarity constrains WIMPs, on the large mass end, to 
be lighter than a few 100 TeV. I think this is a reasonable range to keep in mind, if one 
is wed to the notion of a weakly interacting dark matter particle. 

There exist, however, a number of theoretical prejudices that have populated this 
field for a long time. A historically interesting one has it that WIMPs must be heav- 
ier than 40 GeV. This prejudice somehow even managed to distort how certain dark 
matter search experiments were optimized! It is worthwhile then to see where this 
prejudice comes from. 

The first tenet of the "WIMPs must be heavier than 40 GeV" prejudice is that WIMPs 
are supersymmetric neutralinos. Browsing papers on dark matter (especially on astro- 
ph) the confusion between WIMPs and neutralinos is not unheard of. The second 
tenet is that there exists one universal soft supersymmetry breaking mass for all three 
gauginos at the GUT scale. If Mi(Mgut) = ^2(^gut), where Mi is the soft super- 
symmetry breaking scale associated with the U(l)y gaugino and M2 that of the SU(2) 
gaugino, renormalization group evolution (and the assumption of a "desert" between 
the electroweak and the GUT scale) implies that Mi (Mew) =- 0.4 x M2(Mew)- Now, 
LEP2 constrains the chargino mass to be above about half its center of mass energy, 
or "^vi ~ 100 GeV. But one of the charginos (the wino-like, in SUSY slang) has a 
mass very close to M2(Mew)) therefore implying that the lightest, bino-like neutralino 
m^o ~ Ml (Mew) > 0.4 x 100 GeV = 40 GeV. Amazing. Of course, GUT-scale univer- 
sality, renormalization group evolution etc. are all model-dependent ingredient, not 
to mention the assumption that the dark matter is a neutralino... 

What do we know about the annihilation final state? Presque rien, almost nothing. 
If the dark matter particle is a Majorana fermion, then the pair annihilation into a 
fermion-antifermion final state is "p-wave suppressed": xx ~^ ff requires a helicity 
flip, and thus the matrix element squared is proportional to the square of the fermion 
mass, |Mp oc m'j (in just the same way as for charged pion decay - a Majorana pair 
in a / = wave is in a CP-odd state). As a result, pair-annihilation of Majorana dark 
matter into light fermions is highly suppressed. If < rritop and if the annihilation 
channel xx ~^ bosons (such as W~^W~ or hh) is suppressed, the dominant annihilation 
final states are bb and r+r . This explains the otherwise surprising popularity of these 
two final states in the literature on dark matter indirect detection^. Note that besides 
TJib 3> nir, the bottom quark final state wins by an additional factor 3 from color. There 
exist, however, circumstances where the rr final state can be boosted, for example 
with a light scalar tau in supersymmetry (in the so-called stau coannihilation region) . 

In UED the situation is entirely different: the particle that is usually the stable 
lightest Kaluza-Klein excitation is the n = 1 mode of the hyper-charge gauge boson, 
or The matrix element squared for pair annihilation into a fermion-antifermion 

^If you are a graduate student at UCSC who happened to work on dark matter, and I sit on your 
thesis committee, I will ask you the reason why bb is such a popular final state during your thesis 
defense. Be ready. 
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pair is proportional to the fourth power of the fermion's hyper-charge, |Af p oc \Yf\^, 
thus up-type quarks (y„^ = 4/3) and charged leptons (Y^^ = 2) are the preferred 
annihilation modes. 

If the dark matter lives in an SU(2) multiplet (for example, higgsinos and winos in 
supersymmetry) everything is fixed by gauge interactions. For wino-like dark matter, 
the preferred final state is W^W^. The lightest neutralino is quasi-degenerate with 
the lightest chargino, with mass splittings on the order of a fraction of a GeV, and 
both lie at a scale close to M2, the corresponding soft supersymmetry breaking mass. 
Coannihilation plays obviously a very significant role, and coannihilation is, here, of 
the "symbiotic" type (charginos pair-annihilate quite efficiently) . The resulting pair- 
annihilation cross section is, approximately [21] 

and the resulting thermal relic abundance is 

^ \2.2TeVj ' 

implying that thermal winos must weigh about 2.2 TeV. Higgsinos, instead, come in a 
set of two neutralinos and a chargino with again small mass splittings and at a mass 
scale around /i. The dominant annihilation final states are W~^W^ and ZZ, and the 
pair annihilation and thermal relic densities are given by 

{av) ~ (21 + 3 tan^ Ow + H tan^ Ow) 

512 TT fi^ ^ ' 

and 

indicating that thermal higgsinos like to weigh about a TeV. 
Indirect Detection: Warm-up lap 

Time to charge ahead on indirect detection. The name of the game is to get "enough" 
number counts, in symbols: 

where: 

• 0^ indicates the relevant dark matter-induced event rate, such as the flux of a 
certain type of Standard Model particle, and has units cm~^ ■ s"^; 
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• Acs is an effective area: good to have in mind some numbers here. For exam- 
ple, in the business of gamma-ray telescopes, the Fermi Large Area Telescope 
(LAT) has an effective area of about 1 m^, while the top-of-the-line atmospheric 
Cherenkov (ground based) telescopes, such as H.E.S.S. or MAGIC or VERITAS 
have effective areas on the order of 10^ m^; the relevant numbers for the two 
key antimatter satellites are ~ 0.01 m^ for Pamela and ~ 0.1 m^ for AMS-02; 
finally, if you're asked to quote a number for high-energy neutrino telescopes, 
mention IceCube and mumble 1 km^. 

• Tcxp indicates the relevant "exposure time": for satellites this is on the order of a 
year, which as you know is exactly tt x 10'' s; for typical ground-based telescopes 
you can perhaps count on about lOOh, or about 10^ s, while balloon experiments 
have a typical exposure time of the order of a week, or about 10® s. 

To detect a signal we need to fulfill two basic conditions: 

(i) have some signal events, i.e. 0^ ■ Aes ■ Te^p ^ 1 

(ii) have enough signal-to-noise, for example requesting iVsignai > A/^background- 
I like to classify astrophysical probes of dark matter into three categories: 

1. Very indirect: this category includes effects induced by dark matter on astro- 
physical objects or on cosmological observations; 

2. Indirect: I include in this category probes that don't "trace back" to the anni- 
hilation event, as their trajectories are bent as the particles propagate: charged 
cosmic rays; 

3. Not-so-indirect: neutrinos and gamma rays - we will discuss these guys in the 
next lecture. 

1. Effects on Astrophysical Objects: folks have thought about an amazing variety 
of possibilities, including: 

• Solar Physics (dark matter can affect the Sun's core temperature, the sound 
speed inside the Sun,...) 

• Neutron Star Capture, possibly leading to the formation of black holes (notably 
e.g. in the context of asymmetric dark matter, see e.g. [22]) 

• Supernova and Star cooling (see the excellent book by Georg Rafelt [23]) 

• Protostars (e.g. WIMP-fueled population-Ill stars, available also in Swedish 
[24]) 

• Planets warming 
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Given the relevance of global warming to the general public (and to funding agencies), 
let's make an estimate of this latter effect. The "capture probability" for WIMPs is 
roughly 

'^nucleons ' ^x~^ ' -^planet ^ 10 , 

where for the right-hand side I've used nnucieons ~ Na/chv^, the current upper limit on 
the spin-dependent WIMP-nucleon cross section ct^-n ^ 10^*^ cm^ and the radius of 
Uranus, i? ~ 3 x 10^ cm - the choice of Uranus is motivated by an anomalous heat 
observed in the planet, of about 10^"^ W. Now, the power produced by dark matter 
assuming that all of the dark matter mass is converted to heat is 

W ~ (capture probability) ■ vri^pianet ■ Pdm ■ f dm < 10^^ W 

which tells us that we fall short by a couple orders of magnitude of explaining Uranus' 
anomalous heat. Too bad. 

Exercise #5: Estimate the heat produced by dark matter annihi- 
lation in the Earth and compare with the accuracy of geothermal 
models (see also the much, much more refined discussion in [25]); 
how large should the dark matter-nucleon scattering cross section to 
cause global warming concerns? 



1., cnt'd: Effects on Cosmology: lots of work here, spanning effects on Big Bang 
Nucleosynthesis, on the cosmic microwave background, on reionization, on structure 
formation and many more. I don't even have time to give you a laundry list of all this! 
Go browse the arXiv and have fun! 

2. Charged Cosmic Rays: here, the dark matter "source term" is unfortunately 
tangled with effects of propagation and energy losses of charged cosmic rays on their 
way to our human detectors. 

Do we expect enough cosmic rays from dark matter annihilation or decay to detect 
a signal over the background? The ballpark energy density of cosmic rays in the Milky 
Way is 

ecR ~ J- — ^• 
cni'^ 

Let's estimate the energy density in cosmic rays dumped by dark matter annihilation 
in the Galaxy: 

cdm ~ rriy. ■ {av) ■ n^M ■ ^mw, 
with mx ~ 100 GeV, pdm ~ 0.3 GeV/cm^, {av) ~ 3 x IQ-^e cmVs, and the Milky Way 
age Tmw ~ 10 X 10^ yr, I get 



-2 eV 

eoM ~ 10 



cm^ 



Exercise #6: Improve on the estimate above using a Navarro-Frenk- 
White dark matter density profile and integrating over an appropriate 
cosmic-ray "diffusion region", e.g. a cylindrical slab of half-height 1 
kpc and radius 20 kpc. 
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DM DM ^ atMoM = 1 TeV DM DM ^ W^W atMoM = 1 TeV 




Figure 4: The differential photon (red Unes), neutrino (black lines), (green lines), 
p (blue lines) yield from dark matter pair-annihilation into a qq pair (left) and W~^W~ 
(right). From Ref. [26]. 



Exercise #7: Same as Ex.#6, but for a decaying dark matter parti- 
cle, find tDuiT) where r is the dark matter lifetime. Do you expect to 
get interesting limits on r from this calculation? If yes, please men- 
tion me and these lecture notes in the acknowledgements of your 
forthcoming paper. 



The estimate above indicates that the contribution of annihilating dark matter to 
cosmic rays is, at best, subdominant to the observed cosmic ray energy density, but that 
it could be an 0{1 %) effect. In fact, models of Galactic cosmic rays decently match 
observation, so this is in some sense good news for dark matter model building! As a 
result, it is key in this business to target under-abundant species, namely either heavy 
nuclei or antimatter (for example positrons (e+), antiprotons (p), antideuterons D,...). 
Unfortunately, it is quite hard to produce heavy nuclei from dark matter annihilation 
(that results, in its hadronic part, in a couple of high-energy jets only). Antimatter, 
on the other hand, is promising; typical dark matter models (exceptions are certain 
flavors of asymmetric dark matter) are democratic in producing as much matter as 
antimatter in the annihilation or decay final products. 

Figure 4 illustrates the final yield of several particle species resulting from xx ~^ QQ 
(left) and from xX W^W^ (I took these two nice figures from Ref. [26]). The 
red lines indicate photons, the black lines neutrinos, while the green and blue lines 
indicate e"^ and p, respectively. All of these particle species primarily originate from 
the hadronization and cascade decays of jets initiated by the final state q and q, or 
directly from the prompt decay modes of the W (notice the green and black lines 
getting "horizontal" at x = 1, where x is the particles' kinetic energy normalized by 
the dark matter mass). 
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a. 




kinetic energy [GeV] Energy (GeV) 



Figure 5: The cosmic-ray antiproton to proton ratio (left, from Ref. [27]) and the 
positron fraction (right, from Ref. [28]) as measured by the Pamela experiment. 



In cosmic rays, antimatter is primarily produced by spallation processes, such as 

p+p^p+p+p+p 

where one of the protons in the initial state is a high-energy particle, and the second 
one is typically an nucleus in the interstellar medium gas, and baryon number 
conservation forces you to produce at least four nucleons in the final state. The process 
has a relatively large threshold (if you need a special relativity refresher carry out the 
two-lines calculation), Ep >7 GeV. Now, the spectrum of cosmic rays observed in the 
Galaxy falls steeply with energy. 



diVcosmic— ray protons 7-1— 2 7 



dE 



so compared to the maximal flux of cosmic-ray protons, observed at E ~ 0.1 GeV, 
antiprotons will be under abundant, at 0.1 GeV, by about a factor 



/n 1 \ 2.7 
p I 0.1 



p \7.5 



10 



This is in fact in remarkable agreement with what is observed, see Fig. 5, left, from 
Ref. [27]. 

There are therefore two effects that make antiprotons an interesting probe of dark 
matter (that, as fig. 4 shows, tends to produce low-energy antinucleons) : on the one 
hand there are few "beam" particles to produce cosmic-ray antiprotons, since the 
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cosmic-ray proton spectrum falls steeply, and on the other hand the typical kinetic 
energy inherited by the final state antiproton will be on the same order as the thresh- 
old for the process. Indeed, Fig. 5, left, shows that the p/p ratio peaks right around 
10 GeV, a much higher energy than the typical anti nucleon produced by dark matter. 
These two effects are even more drastic for anti-deuterons (i.e., bound states of p and 
n), for which the key astrophysical background comes from the reaction 

D : p + p p + p + p + p + n + n 

that has a threshold of about 17 GeV. In addition, D have a hard time loosing energy 
by elastic scattering (tertiary population) since the deuteron binding energy is very 
low, and when hit D tend to disintegrate rather than lose energy! There is a smart 
idea out there (the proposed satellite is called GAPS [29]) to target specifically low- 
energy antideteurons and to detect them via the peculiar de-excitation X-rays that an 
atom capturing a D would produce. 



How to deal with charged cosmic rays 

How do we model cosmic-ray transport? The most successful framework is provided 
by the so-called diffusion models (adequate for cosmic-ray energies Eqk < 10^^ eV). 
Let us indicate the differential (in energy) number density of cosmic rays with 

— =ij{x,E,t). 

The master equation of cosmic-ray diffusion models looks something like this: 

-ij = D{E)A^ + —{b{E)^) + Q{x,E,t), (11) 

The first term on the right-hand side describes diffusion, the second one energy losses, 
and the third includes all possible sources. As always in life and in science, it is 
possible (and easy!) to add complications - an incomplete list of popular ones and of 
the associated recipes is: 

• Cosmic-ray convection; recipe: add: ^{vc-ip). 

• Diffusive re-acceleration; recipe: add: ^ Dpp^^ijj. 

• Fragmentation and decays; recipe: add: — 

When dealing with partial differential equations, we all learned in kindergarten that 
it is crucial to define boundary conditions. A popular choice is free-escape at the 
boundaries of a "diffusive region", whose geometry, for obvious reasons, is typically 
chosen to be a cylindrical slab, with 

X 10 kpc, 



— 26 — 



S. Profumo 



Astrophysical Probes of Dark Matter 



h ~ X 1 kpc. 

These numbers (very) approximately reflect the distribution of gas and stars in our 
own Milky Way. 

The diffusion coefficient (that in certain models can depend also on position - it 
more than likely does in reality!) has a dependence on energy (a remnant of the fact 
that the Larmor radius scales with the particle's momentum!) that can be schemati- 
cally cast as 

— , ~ GeV, Do ~ few x 10^^ , 6 ~ 0.7. 

EqJ s 

The parameters entering cosmic ray diffusion are tuned self-consistently to reproduce 
key observational data, such as stable pure secondary to primary ratios as a function 
of energy (classic example: boron to carbon, B/C) or unstable secondary to primary 
ratios, such as ^°Be/^Be. For example, this latter ratio constrains quite severely the 
height of the diffusion region. 

What are the relevant time-scales for the diffusion equations? Two key quantities 
are the diffusion and the energy loss time scales: 

TdiS ~ — ■ 



Do b{E) 

where R is the linear size of the diffusion region, or the relevant time/distance scale 
for which we want to calculate the typical associated diffusion length (for example, to 
infer which diffusion length corresponds to the energy loss time scale, we would plug 
in R ^ c/rioss)- The steady-state diffusion equation (11) can then be re-written as 

= -^ + Q, 

TdiS Tloss 

implying that 

r^Q- min[rdiff, rioss]- (12) 

Let's see if this makes sense and consider cosmic ray protons and primary and sec- 
ondary electrons and positrons: 

• If the primary sources of cosmic-ray protons are supernova remnants, and if the 
injected particles are accelerated via a Fermi mechanism, we expect 



Q ^ E 



-2 



Energy losses for protons in the GeV-TeV range are relatively inefficient, and 
typically Tdig ^ rioss, therefore Eq. (12) would predict 

ip ^ E-'^ ■ E~^ ~ ^-2.7 
which is in great agreement with observation! 
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For primary electrons, let us suppose that again Q ~ E~'^ - for example because 
the acceleration site is the same as for cosmic-ray protons (not such an unrea- 
sonable assumption...). At high energy (Ee GeV) the dominant energy loss 
mechanisms are inverse-Compton scattering (i.e. the process of a high-energy 
electron up-scattering an ambient photon - the inverse of the classic Compton 
scattering where a high-energy photon up-scatters an electron at rest!) and syn- 
chrotron. Both have the same dependence on energy, oc E"^, and the resulting 
energy loss term reads 



where, in units of 10"^^ GeV/s, the constants 

6?c ^ 0.76, ^ 0.025 

and corresponds to the background radiation energy density and B to the 
ambient magnetic field. Depending on the size and geometry of the diffusion 
region, Eq. (12) predicts a break between a low-energy regime where 

''/'primary, low-cncrgy ~ Q ' Tdiff ^ E ^ ■ E ^ ^ E ^''^ 

and a high-energy regime where 

E 2 

^primary, high— energy ~ Q ' ^loss ~ Ej ' ~ 

The general prediction is thus of a broken power-law with a break corresponding 
to Tioss ~ Tdiff- This indeed matches observation again! (both directly and indi- 
rectly, i.e. from measurements of the secondary radiation produced by cosmic 
ray e^). 

For secondary electrons and positrons, produced e.g. by the decay of charged 
pions produced by cosmic-ray proton collisions with protons in the interstellar 
medium, the source term corresponds to the Qp ~ E^'^-'^ spectrum found above. 
The spectrum after diffusion and energy losses will then follow the same 
fate as that of primary particles discussed above: a broken power-law, with a 
hard low-energy spectrum ipsecondary, low-energy ~ E-^-'^ and a softer high-energy 
tail due to energy losses, ipsecondary, high-energy ~ E'^'^ The key point is that, 
independently of the value of 6 (that, remember, tunes the energy dependence 
of the diffusion coefficient) and of the primary injection spectrum and of energy, 
the ratio of secondary to primary species is 

The prediction of a declining secondary-to-primary ratio was recently found to 
be at odds with the observed local positron fraction (see the right panel of figure 
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5, from Ref. [28]), a fact that spurred much speculation about the nature of 
the additional positrons responsible for the upturn in the secondary-to-primary 
ratio. 

There are a couple of special limits in which one can get a simple solution to 
Eq. (1 1) that are worth remembering because they apply to certain interesting physical 
situations: 

1. No diffusion: this case corresponds to physical situations where the energy loss 
time-scale is much shorter than the diffusion time-scale: the cosmic rays effec- 
tively loose their energy before diffusing. In this case, the asymptotic, steady- 
state (i.e. after time transients) solution to the diffusion equation (11) is 



i;(x,E)^-^J dE'Qix,E'). 



There are numerous circumstances where this is a relevant approximation: for 
example, when the system is very large, with its physical size much larger than 
the typical diffusion length associated with the energy-loss time (for example 
clusters of galaxies), or when the system is such that the energy loss term is very 
large (for example, a medium with a very dense radiation fields off of which 
can efficiently inverse-Compton scatter, or with large magnetic field inducing 
large synchrotron radiation energy losses) . 

2. Burst-like injection from a point source at time t in the past: in this case, the 
relevant (spherically symmetric) solution, neglecting energy losses, is 

^ oc Q ■ exp (j-^ j , (13) 

where r is the source distance, 

rdiff - ^/D{E) ■ t. 

This second case, a burst-like injection, is especially important in connection with 
Galactic pulsars as sources of high-energy e^, a potential explanation to the anomalous 
rising positron fraction found by Pamela [28] and recently confirmed by the Fermi-LAT 
[30], as pointed out by many (including Yours Truly, who once again is not shy to self- 
promote his own papers, see e.g. [31]). 

What are the requirements on the age and distance of a pulsar that contributes to 
the Pamela positron anomaly? It is easy to give general arguments: first, the pulsar 
age must be much shorter than the energy loss time scale for energies as large as 
about Ee ~ 100, in order to have at least some energetic around! This implies the 
condition 

Tpsr < Tioss = for E = 100 GeV, n,,, ~ iq-1^^iqq2 ^ ~ ^O'" ^ ~ ^ Myr. 
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Now, to avoid the exponential suppression of Eq. (13) we must have 



D{E) ■ Tpsr > distance ^ distance < (3 x lO^^-lOO^'^-lO^^)^/^ cm ~ 10^^ cm ~ 3 kpc. 



So our candidate pulsar is younger than about a mega-year and closer than a few 
kilo-parsec. One would also like the pulsar to have enough power injected in electron- 
positron pairs, but this condition, for such a nearby object, is usually fulfilled. Inter- 
estingly, several pulsar candidates exist within the desired age and distance, including 
possibly a handful of the newly discovered radio-quiet gamma-ray pulsars detected by 
Fermi-LAT (see e.g. Ref. [32]). 

As some of you might be aware of, the dark matter annihilation explanation to 
the Pamela positron fraction anomaly gathered quite a bit of attention (in fact on the 
order of 10^ publications entertain this possibility!). Dark matter as a source of the 
observed excess high-energy positrons faces various issues, including the following: 

• there is no evidence for an associated antiproton excess, thus the dark mat- 
ter must preferentially pair-annihilate into non-hadronic final states (it must be 
"leptophilic"); 

• diffuse secondary radiation from internal bremsstrahlung and inverse-Compton 
is not observed; 

• the needed pair-annihilation rate. 



is very large for thermal production, and generically leads to unseen gamma-ray 
or radio emission; 

• a XIII century monk pointed out that "entia non sunt multipUcanda prseter necessi- 
taterri", and the pulsar explanation works just fine to explain the excess positrons. 

Despite these difficulties, theorists from all over the world (including myself) have 
proposed models that circumvent all difficulties and show proof that Pamela might 
have perhaps detected the first non-gravitational signs of dark matter, providing more 
and more empirical evidence in favor of the so-called "Redman theorem" [33] : 

Any competent theoretician can fit any given theory to any given set of facts 

What's next? Well, myself and many others are quite eagerly awaiting results 
from the AMS-02 payload, successfully deployed and operational on the International 
Space Station since May 2011. An independent measurement of the positron fraction 
over an extended energy range (especially in the critical high-energy end, where a 
cut-off in the positron fraction might indicate a new physics origin!) and with much 
larger statistics, measurements of various cosmic ray species which will be key to a 
better understanding of cosmic ray propagation in the Galaxy, and possibly additional 
information on e.g. anisotropies of the positron arrival direction don't quite keep me 
awake as much as my newborn son (actually, not even nearly), but still. . . 
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Lecture 4: not-so-Indirect Detection: Neutrinos and Gamma 
Rays 

The tiny neutral ones 

Detecting neutrinos (from an Italian made-up word that indicates the "tiny neutral 
one", with an English made-up plural form) is hard. In fact, despite building km^ 
size detectors, only two astrophysical neutrino sources have been observed so far: the 
Sun and Supernova 1987A! The flip side of the coin is that astrophysical backgrounds 
are evidently quite low (albeit of course cosmic rays produce copious "atmospheric" 
neutrinos as they hit the atmosphere...) if one is to use neutrinos to search for dark 
matter. The key idea is that dark matter particles can accrue in celestial bodies until 
large enough densities start fueling a steady rate of annihilation yielding high-energy 
neutrinos. Neutrinos are pretty much the only thing produced by dark matter annihi- 
lation that can escape the core of a celestial body without losing much energy at all, 
and get all the way out to our km^ size detectors. The best bets are the Sun and the 
Earth, with the former, turns out, much better than (although somewhat complemen- 
tary to) the latter. Let us now make a few estimates for this process, for the case of 
dark matter capture and annihilation in the Sun. 

The dark matter capture rate in the Sun is, roughly 



with the dark matter flux 



"x 

the ratio Mq/uip estimating the number of target nucleons in the Sun, and the dark 
matter-nucleon interaction cross section a^-p being bound by current experimental 
limits: 

spin dependent < -i rv-39 2 
X~P ~ ' 

spin independent ^ -i n— 44 2 

(Tt^ — p J-^ cm . 

Plugging in the relevant numbers, I find 

^ s V0.3 GeV/cm3 y ' V300 km/s J ' \ )' VlQ-^s cm^/ 

We are interested in the number of dark matter particles in Sun: let's call this number 
N and write down a differential equation that describes the time evolution of N{t): 

^ ^ - A'^lNit)? - E^N(t), 
at 

There are various elements I introduced in the otherwise self-explanatory equation 
above: 
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describes the "evaporation" of dark matter particles, something that happens 
if the particles have a (thermal) velocity comparable with the celestial body's 
escape velocity. Let's quickly estimate this effect. For the Sun 

o km o 
i;® ~ 1156 3 X 10"^ c 

esc g 

while the Sun's core temperature (the dark matter particles sink to the center 
after multiple scattering inside the Sun) is 

Tg^^^m' K^lkeV^m^-vl. 

This gives, for the typical dark matter thermal velocities in the core of the Sun 



1 keV 



1/2 



c • y—j Z €c ^rn^< 0.1 GeV. 

Bottom line: for dark matter particles in the "preferred" WIMP mass range we 
can safely neglect evaporation. 



Exercise #8: Estimate evaporation in the case of the Earth: 
what is the relevant dark matter particle mass range for which 
evaporation matters in the Earth? 



The annihilation rate 

where Ves is an effective volume which depends on where WIMPs live inside the 
Sun; let us use the following (reasonable) guess for the density profile of the 
sunk WIMPs in the Sun: 



n(r) = no exp 



;;0grav(''") 



J'0 



We can choose to estimate by identifying an effective radius it!efr correspond- 
ing to the condition 

0grav(^eff) _ . _ GjvP® f i^eff^X 



^ i?efr ~ lO'' cm ( ) ^ Kff ~ 10'" cm=^ (- 

V 100 GeV/ v: 



"x 



-V 



.100 GeV/ ' V 100 GeV/ 

Remember that the Sun's radius is approximately i?® ~ 7 x 10^° cm, so this 
radius is smaller than the Sun's radius for reasonably light WIMPs. 
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Neglecting evaporation, even I can solve the differential equation above, and calculate 
the quantity we are really interested in, the annihilation rate in the Sun: 



r^ = ^A®[iv(t®)]2 



with 



~ 4.5 Byr ~ 10'^ s 



tanh(\/C®A® i®) 



17 



the Sun's age (not to be confused with the Sun's core temperature T®!). One thing we 
learn from the solution above is that equilibrium between capture and annihilation is 
reached if 

Do we expect equilibrium or not, for nominal WIMP parameters? Yes, we do! Let's 
plug in the numbers and convince ourselves of this fact: first, let's find the required 
annihilation rate for equilibrium 



1 



x-p 



10-39 



4® ^ 

eq ^ (^0^2 C® 1034 . 1023 g 



10 



Now, for vanilla WIMP dark matter 

A® 



3 X 10"^^ s"^ 



(av) 



3 X 10-26 cmVs 
so equilibrium is reached for small as about 10 cra^ 



Exercise #9: Re-do this calculation for the case of the Earth and find 
the critical dark matter-nucleon scattering cross section for equilib- 
rium; note that the relevant scattering cross section in the Earth is 
spin-independent (as the Earth is mostly made of spin-0 Iron nuclei) : 
do you then expect the equilibrium condition to hold for the flux of 
neutrinos from the center of the Earth? 



If equilibrium is achieved, then 



c® 



and we don't care about the pair annihilation cross section (a unique case in the 
business of indirect dark matter detection!), while we only care about the cross section 
for dark matter capture. The resulting flux of neutrinos of flavor / will then be 



dA^, 



C® 



dE 87r(L>®)2 \dE, 



dN, 
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where the last factor with the subscript "inj" is the "injection" spectrum of neutrinos 
per annihilation. Effects that complicate this discussion include neutrino oscillation, 
absorption of neutrinos in the Sun, and many others that a few smart people out there 
have already kindly worked out for you. 

The final step is to count the number of events we expect at IceCube or at any other 
mega-neutrino-detector (these detectors are fundamentally arrays of photomultipliers 
reading Cherenkov light from muons produced by z/^ charged-current interactions) : 

iVevents = J dE,^ J dy (^A,s " ^ " ^{E,,, y) " (^m(^^J)) , 

where y indicates the energy fraction transferred to the /z in the charged current 
interaction, da/dy is the relevant cross section for charged current interactions, and 
the last factor i?^ indicates the muon range in the relevant material the detector lives 
in (for example, Antarctic ice for IceCube). 

The most promising dark matter pair annihilation final states in this business are 
those producing a "hard" spectrum of muon neutrinos, i.e. energetic neutrinos. These 
by all means include W^W^ and ZZ pairs, that dump out prompt muon neutrinos 
from the leptonic decay modes of the gauge bosons; luckily, for example in supersym- 
metry these are exactly the preferred final states for wino- and higgsino-like lightest 
neutralinos. 

The typical flux sensitivity threshold we want to hit to get an interesting signal is 
about hundreds of muons per km-squared per year, and the typical energy thresholds 
are 100 GeV for IceCube, which is improved down to 10 GeV for DeepCore and that 
could go down to the order of a GeV for the further thickly instrumented portion of 
the detector to be named PINGU. 

Light from Dark Matter 

There are two key ways to get light out of dark matter: 

(i) Prompt photons from the annihilation or decay event, and 

(ii) Secondary photons from radiative processes associated with stable, charged par- 
ticles produced by the dark matter annihilation or decay event (in practice, the 
most important ones are electrons and positrons) 

Prompt photons are produced either by the two-photon decay of neutral pions 
7]-o _i. dumped by the hadronization chain of strongly interacting annihilation prod- 
ucts, or by internal bremsstrahlung off of charged particles in the intermediate or final 
state; this second contribution is typically "harder", i.e. more energetic, than the first 
one. Gamma rays from neutral pion decay have the nice spectral feature that I ask 
you to derive in the next exercise^. 

remember this problem well, as it was asked to me during my PhD entrance exam by my advisor- 
to-be! 
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Exercise #10: Show that, independent of the 7r° spectrum, the dif- 
ferential spectrum of gamma rays resulting from 7r° — > 77, dN'!^'^ /dE^ 
is symmetric around E^^ = on a log scale in energy. 



Secondary photons originate as the counterpart of the key energy loss processes 
for electrons and positrons we discussed in the previous lecture: inverse-Compton and 
synchrotron. To qualitatively understand the features of inverse-Compton emission, 
it is useful to commit to memory the formula for the average energy {E'q) of the up- 
scattered photon (with an original initial energy E^) as a function of the Lorentz factor 
7e = Ee/nie of thc impinging high-energy electron: 

The relevant numbers for Eq are as follows: 

CMB : ~ 2 X 10"^ eV 
starlight : i^o ~ 1 eV 
dust : Eo ~ 0.01 eV 

so for a typical electron-positron injection energy from dark matter 

Ee ^ ^ 7e ~ 2 X lOM ^ ) 

10 ' V 100 GeV/ 

and 

Inverse-Compton emission from dark matter therefore produces hard X-ray photons 
in the hundreds of keV range. This is great news, as a brand new NASA telescope, 
NuSTAR, is looking at the sky exactly in that energy range [34] ! The inverse-Compton 
light from starlight and dust falls, instead, in the low-energy gamma-ray regime. 
In the monochromatic approximation, synchrotron emission peaks at 

,, / F \ / R \ 



^sync ^ 



MHz V GeV 



2- ^ 



) UgJ 



and the synchrotron power scales like 5^. Dark matter annihilation thus produces a 
rich, multi-wavelength emission spectrum that goes well beyond the gamma-ray band. 
An example of the spectrum expected e.g. from the nearby Coma cluster of galaxies is 
shown in Fig. 6, left, from Ref. [35]. Note that the various secondary emission peaks 
appear exactly where the formulae above would predict them to be! 

While secondary emission is always present, it involves the additional steps of 
accounting for the diffusion and energy losses of the produced by dark matter 
annihilation. Prompt gamma-ray emission, on the other hand, is simpler, and it only 
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Figure 6: Left: The multi-wavelength emission spectrum from the pair-annihilation 
of a dark matter particle with = 40 GeV in the Coma cluster of galaxies, from 
Ref. [35]. Right: the pair annihilation cross section into two photons for MSSM neu- 
tralinos, from Ref. [37]. 



involves identifying a dark matter structure and a particle dark matter model; we 
will thus here make a few estimates for this prompt emission only, which for nominal 
WIMPs produces photons in the gamma-ray energy range. Also, for definiteness we 
will talk about dark matter annihilation - dark matter decay is even simpler! 

What are the optimal targets and the expected detection rates for gamma ray 
searches for dark matter? The flux of photons produced by dark matter annihilation 
from a given direction ip in the sky and from within a solid angle is 

where the last factor is, in fact, a sum of the prompt gamma-ray spectrum from every 
annihilation final state /: 

dE^ ~ ^ dE^' 

and where the term in curly brackets is often referred to as the "J-factor", and is a 
function of the solid angle Af2 and of course of the direction in the sky, J = J(Af2, 
and carries units of GeV^/cm^. The solid angle A^l should be optimized for a given 
gamma-ray detector and for a given target, field of view and angular resolution to 
maximize (typically) the signal to noise. It turns out that the relevant solid angles 
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correspond to an angular extent of about one degree, or Ail ~ 10 sr for the Fermi- 
lAT at an energy of about a GeV, down to an angular extent of 0.1 degrees, or Afl ~ 
10~^ sr for ACT, or for Fermi in the high-energy regime. 

Let me give you a "laundry list" of potential interesting targets to search for a 
gamma-ray signal from dark matter; for most of these targets the "J factor" is approx- 
imately the same for a solid angle corresponding to 1 deg or 0.1 deg: 

1. Dwarf Spheroidal Galaxies 

• Draco, J ~ lO^^ GeVVcm^ ± a factor 1.5; 

• Ursa Minor, J ~ 10^^ GeVVcm^ ± a factor 1.5; 

• Segue, J ~ 10^0 GeV^/cm^, ± a factor 3 

2. Local Milky-Way-like galaxies 

• M31, J ~ 10^° GeVVcm^ 

3. Local clusters of galaxies 

• Fornax, J ~ 10^^ GeVVcm^ 

• Coma, J ~ 10^^ GeVVcm^ 

• Bullet, J - 10^^ GeVVcm^ 

4. Galactic center 

• 0.1°: J ~ 10^2 1025 GeVVcm^ 

• 1°: J ~ 10^2 . . . 10^4 GeVVcm^ 

To have a detection, we need to have enough photon counts, possibly a lot: 

J E~i range 

The following table gives a rule of thumb for the relevant energy ranges, effective 
areas and observing time for current and future gamma-ray observatories: 





Fermi-LAT 


H.E.S.S. 


CTA 


E^ range 

T 

obs 


0.1 to 300 GeV 

~ 1 

~ 10^ s 


0.1 to 10 TeV 

~ 10^ 
~ 10^ s 


10 GeV to 10 TeV 

~ 10^ m2 
~ 10*^ s 



It is instructive to calculate the minimal J factor needed to get at least some 
gamma-ray signal from dark matter. Consider for example Fermi-LAT: over the LAT 
energy range, typically 
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so that 



X 

and 



^ ' M w? ^ cm2 s VGeVVcmS J 



J 



GeVVcmS' 

where we put in the nominal values for the effective area and observing time as in the 
Table above, so, we want 

J > 10^° GeVVcml 

This is a bit bigger than the individual dwarf spheroidal galaxies' J factors quoted 
above. In fact, combining observations of all (non-detected) dwarf galaxies gives one 
of the tightest (in my personal opinion the tightest) constraints to date on the dark 
matter pair-annihilation rate as a function of mass: dwarf galaxies are a virtually 
background free target, with 

J'°* ~ few X 10'° GeVVcm^ 
therefore the resulting limits [36] are 

cm^ /30 GeV 



(cTf) lim ~ 3 X 10 



-26 



A clear-cut signal of dark matter annihilation that could be detected with gamma-ray 
detectors is the close-to-monochromatic line from the direct annihilation 



XX 77- 

Since ~ m^, ~ and the line is almost monochromatic. Dark matter particles, 
being dark, ought better not directly couple to photons, and the naive expectation is 
that the 77 amplitude be loop-suppressed, i.e. 

rsj ^ 

(crv)tot 167r2' 

This naive estimate gives the correct ballpark over a wide range of parameters in the 
MSSM, as fig. 6, right (taken from Ref. [37]), illustrates (but things can get a lot 
worse, and not much better, than the naive estimate!). As you might have heard, 
a monochromatic line at 130 (or perhaps 135) GeV might indeed be present in the 
Fermi LAT data [38, 39], although a solid confirmation (possibly only with the new 
generation LAT analysis software Pass 8) that this is not an instrumental effect, and 
that the signal is indeed as statistically robust as it (frankly very much) looks, has yet 
to come. Stay tuned. Clearly, very exciting days lie ahead! 
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